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The dark matter problem will be solved only when all of the dark matter is ac- 
counted for. Although wimps may be discovered in direct detection experiments 
soon, we will not know what fraction of the dark matter halo they compose until 
we measure their local density. In this talk, I will offer a novel method to determine 
the mass of a wimp from direct detection experiments alone using kinematical con- 
sistency constraints. I will then describe a general method to estimate the local 
density of wimps using both dark matter detection and hadron collider data when 
it becomes available. These results were obtained in collaboration with Gordon 
Kane at the University of Michigan. 



1. Introduction 

The direct detection of wimps in the galactic halo would be an enormous 
triumph of experimental and theoretical particle cosmology, have deep im- 
plications for our understanding of the universe, and would explain (at 
least) some of the dark matter in the universe. However, it is unreasonable 
to assume that the entire dark matter halo is composed of only the wimp 
observed. 

What fraction of the dark matter halo is represented by a particular 
wimp is a question that cannot be answered by direct detection experi- 
ments alone or colliders alone. Even if a stable, weakly interacting massive 
particle is discovered at the LHC, no amount of collider data can measure 
its actual local density. Alternatively, even if many direct detection experi- 
ments unambiguously observe the scattering of wimps from the halo, these 
experiments cannot identify the wimp or determine its couplings — crucial 
to measuring the local density. This point was first raised technically in 
Ref."'^. In this talk, I will present a general way to address the cosmological 
significance of wimps observed in direct detection experiments using data 
and bounds from colliders. For a more detailed discussion, see Ref.^. 
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2. The No-Lose Theorem vs. Our Ability to Win 

Let us imagine that a weakly interacting massive particle x has been un- 
ambiguously observed in direct detection experiments. There is no reason 
to suspect that x is o,ll the dark matter. In fact, it is very easy to find 
models where detectable wimps compose only a fraction of a percent of all 
the dark matter. Using the framework of the DarkSUSY code'', wc ran- 
domly generated some six thousand constrained MSSMs and found that 
for any signal rate, the relic density fluctuates over at least two orders of 
magnitude^*^. It appears as if the detection signal rate and relic density 
are largely uncorrelated: a wimp with a low relic density may be just as 
observable as one making up most of the dark matter. 

The fact that particles with even small relic densities can have large 
detection signals has been noted by many authors (sec, e.g. Rcf.^) and has 
sometimes referred to as the 'no-lose theorem:' experimentalists may not 
lose out on discovering even a very tiny faction of the dark matter halo. 
However, the no-losc theorem also implies that a wimp discovery could eas- 
ily represent a negligible fraction of the dark matter. Therefore, although 
wimps may be discovered in the near future, the dark matter problem will 
not be solved until the density of wimps has been directly determined, and 
all the dark matter is accounted for. 

3. Dgirk Matter Direct Detection Rates 

In general, the wimp-nucleon elastic scattering scattering rate is a function 
of the cross section for scattering, nuclear physics describing the detector, 
and the local velocity profile of the wimp fraction of the dark matter halo. 
If a detector is composed of nuclei labeled by the index j, each with mass 
fraction Cj, then the differential rate of wimp scattering at recoil energy q 
is given by''. 



"Only an upper bound on the relic density, consistent with WMAP, was imposed while 
generating these models. 

''A detailed discussion of equation 1 can be found in most modern reviews of dark matter, 
(see, e.g., Ref.^). 
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Figure 1. The function C,(m'^) where the wimp corresponds to the neutrahno in the 
MSSM specified by ATLAS SUSY point 2. The models and data were generated within 
the framework of the DarkSUSY package. 

where f{v,t) is the halo velocity profile, Fj{q) and Sj^^{q) are nuclear 
form factors, oq = Cp + a„ and ai = Up — Un, and the constant parameters 
fp^n and Up^n describe the coherent and incoherent wimp-nucleon scattering 
cross sections, respectively. 

From the equation above, it is clear that to determine the density p^, 
one must identify the particle, determine its mass, estimate the halo profile, 
and 'know' the interaction parameters from the theory describing x- Each of 
these require enormous efforts of both dark matter and collider experiments. 

4. Determining the Wimp Mass 

Initially, perhaps the most important parameter of x to determine is its 
mass. Not only does determine the particle's kinematics, but it may be 
critical to the identification of the particle'^. 

There are two known ways to determine from direct detection data 
alone. One method, using the annual modulation crossing energy, was de- 
scribed in a dark matter review article by Primack et. al. in 1988^'*. The 
other method has been developed by the author and is described presently. 

If the halo velocity profile and m-^ are known, then direct detection 
data from different detector materials and different energies can be used to 



'^One may hope that accelerators observe a weakly interacting, stable, massive particle 
with the same mass as x- 

Although it seems unlikely to have originated in a review article, we have been unable 
to find any earlier reference. 
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solve for ^/p^fp,n and y/p^ap.n by inverting equation 1. If the halo velocity 
profile can be approximated, only is required to compute y^fp,n and 
y/Pxap,n with sufficient data. We can generally expect to have many more 
measurements than the minimum required to solve the system of equations 
once wimps arc observed. 

Because the interaction parameters are constant, all linearly indepen- 
dent combinations of measurements used to solve for the scaled interaction 
parameters must agree, if the correct mass were used in the derivation. 
This motivates us to define a 'kinematical consistency' function ({m'^), 

C("^x) = a/Px H {^"p(*) ~ + - an{j)f + similar termsj , 

which compares the values of ^/p^jp.n or y/p^ap.„ obtained using different 
independent subsets of the data — indexed by i, j — as a function of used 
to invert the equations. It is obviously necessary that C(™x) ~ ^ when 

To determine the wimp mass, one varies m!^ until C,(rn'^ = 0. We ap- 
plied this test to some six thousand random, constrained MSSMs. For every 
single model tested, the correct mass was determined to near-arbitrary pre- 
cision. Figure 1 illustrates a typical plot of C,{m'^). Notice that it has an 
extremely sharp minimum within a few GeV of the true wimp mass. It 
should be noted, however, that experimental uncertainties and resolutions 
were not considered during these calculations. 

5. Neutralino Dark Matter and a Strict Lower Bound on p,^ 

Before we can compute p^, the interaction parameters /p^„ or ap^n must be 
'known.' These interaction parameters depend on very detailed knowledge 
of the particle physics of x- Let us consider the specific case in which the 
discovered particle is the neutralino. The interaction parameters will then 
depend on many details of the MSSM — these may not be known until well 
after dark matter particles have been discovered. Nevertheless, even with 
extremely limited knowledge of the MSSM, we can estimate the parameters 
using partial data, bounds, and constraints. 

For example, we have found that given bounds on tan/3 and a lower 
bound on the lightest squark mass, m^, there is a strict upper bound for 
the incoherent x-quark scattering parameters. In this case, it can be shown. 



•^This may not be a sufficient condition, however; although, we have found no example 
where Cim'-^) = when rra^ ^ m^. 
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that the magnitude of a„ is strictly bounded by 

a« < - -J^iNl -Nl) + 7^ ^ tan^ 9wNI + ^M, 

777-^ 1 

+ 2 "2^ ^i +-tan0iy|Ar5||Ar^|cos(a^) 

77? 77? 1 

+ ^\Nyir\\N{j Icosfas -a^) ^!—-ta.n0w\Nf,\\Nfv |cos(a^ ) } 

mvKsin/?/ ^^2' v mw sin 0e ' ^" -f^^' ^ -f^^^J 

where a^^^, and are the relative phases between N^^,N^ and N^, 
respectively. This expression has six real unknowns. Notice that by the 
normalization of the neutralino wave function, the parameter space is com- 
pact. Therefore, a„ can be absolutely maximized with respect to all six un- 
knowns. It should be emphasized that this analysis is for the most general 
softly-broken MSSM; no ad hoc supersymmetry breaking scenarios — such 
as mSUGRA — were assumed. 

Because we can determine ^/p^Op^n, the strong upper bounds on ap,„ 
directly translate into strong lower bounds on . To test this idea, we con- 
sidered some six thousand randomly generated, constrained MSSMs. For 
each of these models, upper bounds were calculated for Op,„ assuming 10% 
uncertainty in tan /? and a lower bound on the lowest squark mass of either 
200 GeV or the actual mass of the lightest squark, whichever is less. The 
specific gauge content of the neutralino was taken to be known for each 
model for computational simplicity^. Using the upper bounds for ap,„, we 
obtain a lower bound on the local density p^. 

Figure 2 illustrates the results of using this algorithm. Notice that the 
estimated local density is always strictly less than the true local density. 
Also, for many models the lower bound is not such a poor estimate. 

6. Conclusions 

We have seen that, by itself, a discovery of dark matter particles in our 
galactic halo cannot solve the dark matter problem. However, combined 
with data from colliders to identify a particle and determine its interaction 
parameters, we can generally estimate its local density. 

We presented a robust, model independent method to determine the 
mass of a wimp using direct detector data alone. We have shown explicitly 



If the gauge content of the neutrahno was unknown, the interaction parameters could 
have been maximized with respect to tliesc parameters as described earlier. In general, 
therefore, the upper bounds plotted are more restrictive than they would be in practice. 
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Figure 2. This plot compares the lower bound and estimate of the local denisty com- 
puted using the strong upper bound for ap_„ to the true local density for each model. 

how one can determine or estimate in the case where x is the neutraUno. 

Ahhough the dark matter problem may not be solved immediately when 
wimps are discovered, there are clear and general ways to address their cos- 
mological significance. 
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